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Abstract 



f — ■ We performed simultaneous spectral analysis of 26.6 ksec of Chandra and 102.3 ksec of Suzaku X-ray 

£NJ ', data of the starburst galaxy NGC 3079. The spectra are extracted from four regions: 0.5' (2.25 kpc) circle, 

an inner 0.5'-l' (2.25-4.5 kpc) ring, and an outer l'-2' (4.5-9 kpc) ring from Chandra, and 4' (18 kpc) 
circle from Suzaku, all centered on the nucleus. Fitting with thermal plasma models yields interstellar 
medium (ISM) temperatures of 0.65tp;°4, 0.45to;°g, and 0.24l[];^ keV in the three regions, respectively. 
' The combination of Chandra's high angular resolution and Suzaku's good spectral sensitivity enable us 

to spatially resolve and measure the abundances of the metals O, Ne, Mg, and Fc within the hot ISM. In 
particular, the abundance patterns of O/Fe, Ne/Fe, Mg/Fc, and Si/Fe in the central regions (<4.5 kpc) 
are consistent with the expectations from Supernovae (SN) II synthesis. On the other hand, the pattern in 
the region beyond 4.5 kpc is closer to solar. The central regions are also where copious polycyclic aromatic 
hydrocarbon infrared emission related to recent starburst activity is known to occur. This suggests that 
. we are seeing starburst-related SN II metal enrichment in the hot X-ray-emitting nuclear ISM. The spatial 

extent of the SN II— like abundance patterns is consistent with NGC 3079 being in a relatively-early phase 
of starburst activity. 

Key words: galaxies: individual (NGC 3079), galaxies: interstellar medium, galaxies: abundances, 
| galaxies: starburst, X-rays: galaxies 

o 

' 1. Introduction lar abundance (Loddcrs 2003), respectively (Yamasaki et 

• , al. 2009). These results support the scenario where star- 

Thc hot interstellar medium (ISM) in galaxies hold im- burst activity enriches the outflow through SN II metal 
portant information to reveal the chemical evolution of ejection into intergalactic space. In contrast, NGC 4258 is 
galaxies and clusters of galaxies. O and Mg are synthe- a non-starburst spiral galaxy, whose ISM abundance pat- 
sized by type II supernovae (SN II), while Fe is mainly tern is fully consistent with solar (Konami et al. 2009). In 
. ' produced by type la supernovae (SN la). Therefore, the fact, spiral galaxies at the present cosmic epoch commonly 
ratios of metal abundances provide us clues to investigate show solar abundance ratios unaffected by starburst ac- 
■ the star formation history in galaxies. In addition, star- tivity. 

burst galaxies play a prominent role in the chemical evo- NGC 3079 is a nearby SB(s)c starburst galaxy with a 
lution of the universe since outflows inject metals into the high inclination angle of 84° and a disk position angle of 
intra-cluster medium (ICM) and warm-hot intergalactic 166° (Irwin & Seaquist 1991) determined in the optical, 
medium (WHIM). Therefore, this galaxy gives a clear view along the outflow 

Observations of starburst galaxies in X-rays have re- axis perpendicular to the disk, suitable for investigating 
ccntly been performed to measure the metal abundance chemical enrichment by the starburst. The derived radial 
patterns in the ISM. XMM-Newton detected the metals velocity of the CO emission line imply the existence of an 
O, Ne, Mg, Si, and Fe at < 3 kpc from the nucleus along ultra-high-density core of molecular gas within a radius 
the outflow axis of NGC 253 (Bauer et al. 2007). For of 125 pc, with mass as large as ~ 3 x 10 s M Q (Sofue et 
M 82, the abundance patterns were derived in the outflow al. 2001). In the far-infrared, the luminous and extended 
to 'cap' regions located ~11 kpc to the North (Tsuru et emission is distributed across the entire disk of the galaxy 
al. 2007; Konami et al. 2011). The ratios of O/Fe, Ne/Fe, rather than being concentrated in an central region (Perez 
and Mg/Fe are higher than the solar values by factors Garcia et al. 2000). The observations of polycyclic aro- 
of 2-3, as expected from SN II yields. Furthermore, for matic hydrocarbons suggest the occurrence of active star- 
NGC 4631, the outflow wind and disk regions show pat- formation in the central 4 kpc region (Yamagishi et al. 
terns close to those expected from the SN II yields and so- 2010). The outflow with kiloparsec-scale along the minor 
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Fig. 1. The left and right-hand panels are the X-ray images of NGC 3079 with Chandra ACIS-S3 and with Suzaku XIS1, respectively. 
Background components were not subtracted and vignetting was not corrected for in these images. These images are overlaid with 
optical galactic disk contours from the Digitized Sky Survey (6450A) in green. Left; the ACIS-S3 image in 0.3-5.0 keV is smoothed 
with a ct=1.95" Gaussian profile. The radii of the white circles are 0.5', 1', 2', and 3' moving out from the nucleus. The white 
box shows the ACIS-S3 chip boundary. Right; the XIS 1 image in 0.5—5.0 keV is smoothed with a cr=6.24" Gaussian profile. 
The radius of the white circular source region is 4'. Background spectra were extracted from the full Suzaku XIS field of view 
excluding the source aperture and the two yellow circles centered on field sources. The white box is the Suzaku XIS field of view. 



axis have been observed in the radio (Duric & Seaquist 
1988; Sofue et al. 2001), Ha (Ford et al. 1986; Veilleux et 
al. 1994), and X-ray emission (Dahlem et al. 1998; Pietsch 
et al. 1998). In particular, there is a tight correlation be- 
tween Ha and X-ray emission in the superbubble, but 
the nature of the power source of the superbubble is still 
debated, with both active galactic nucleus and starburst 
sources discussed (Cecil et al. 2002). Pietsch et al. (1998) 
reproduced the spectra of ROSAT X-ray data in the en- 
ergy range of 0.1-2.4 keV by using thermal plasma models 
with a negative radial temperature gradient moving out- 
wards. 

NGC 3079 is also classified as a low-ionization nuclear 
emission-line region or as a Seyfert 2 galaxy in the optical. 
Iyomoto et al. (2001) detected a heavily absorbed contin- 
uum and a strong Fe-K emission line, using BcppoSAX 
data in the energy range 2-100 keV. They estimated that 
the contributions of the AGN and the starburst to the 
energy output of NGC 3079 are of the same order of 
magnitude. Fukazawa et al. (2011) systematically ana- 
lyzed the Fc-K line features of Seyfert galaxies including 
NGC 3079with Suzaku data. They revealed neutral mat- 
ter causing a 6.4 keV line and 7.1 keV edge is attributed 
to the torus, constraining the ionization state and the re- 
lation among lines, column density, and edge depth. 

In this paper, we investigate the spatial abundance pat- 
terns in the hot ISM of NGC 3079. The paper is struc- 
tured as follows. In section 2, we summarize observations 
by Chandra and Suzaku. Chandra's sharp angular reso- 



lution allows spatial resolution of emission from the ISM 
and individual contaminating point sources. Suzaku's ex- 
cellent spectral sensitivity allows high quality modeling of 
the abundance patterns. Section 3 details the data anal- 
ysis and results, and section 4 gives a discussion of these 
results. Finally, in section 5, we present our conclusions. 
We adopted 15.6 Mpc for the distance to NGC 3079 (Sofue 
et al. 1999), where 1' corresponds to 4.5 kpc. Unless noted 
otherwise, we use the solar abundances from Lodders 
(2003), and the quoted errors are for a 90% confidence 
interval for a single interesting parameter. 

2. Observation and Data Reduction 

2. 1 . Chandra 

Chandra observed NGC 3079 (ObsID 2038) in 2001 
March, using the Advanced CCD Imaging Spectrometer 
(ACIS; Garmire et al. 2003). The aimpoint on the default 
location of the S3 chip is (R.A., Dec.) = (10 h 01 m 53f69, 
+55°40'41"1). We reprocessed all the level-1 event data, 
applying the standard data reduction using CIAO ver- 
sion 4.3 and CALDB version 4.3.3. The light curve over 
0.4-10 keV are extracted with CIAO tool "dmextract". 
There were no anomalous events at greater or less than 
3cr around the mean count rate, so the full exposure time 
of 26.6 ksec was retained. 
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Table 1. Best-fit parameters for the integrated spectrum of 
point sources (including the nucleus) of Chandra with apec 
components + PL model.* 

Parameters 



^H PS 

Tps 

line energy 
X 2 /d.o.f. 



(10 20 cm- 2 ) 
(kcV) 



ii0 -35.4 

1 Q2+ - 15 

°- JU -0.03 

46/36 



Spectral analysis was performed with HEAsoft ver- 
sion 6.10. In order to subtract the non- X-ray back- 
ground (NXB), we employed the dark-Earth database 
available through the "xisnxbgen" ftools task (Tawa et 
al. 2008). We generated two different Ancillary Response 
Files (ARFs) for the spectrum of each region, which as- 
sumed the source image as uniform sky emission and the 
observed XIS1 image by the "xissimarfgen" ftools task 
(Ishisaki et al. 2007). We also included the effect of con- 
tamination on the optical blocking filters of the XIS in the 
ARFs. 
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Fig. 2. Background-subtracted Chandra ACIS-S3 spectrum 
(black data) of accumulated emission for point sources, 
shown without removal of instrumental response. The black 
histogram shows the best- fit model from table 1. The 
green line shows the absorbed power-law (PLps) component, 
and the purple line is the Fe-line component (Gaussian). 



2.2. Suzaku 

Suzaku performed an observation of NGC 3079 in 2008 
May, with the target centered on the HXD aimpoint 
(ObsID 803039020). The averaged pointing direction is 
centered on (R.A., Dec.) = (10 h 01 m 59f33, +55°36'41."8). 
The XIS consists of two front-illuminated CCD cameras 
(FI: XIS0 and XIS3), and one which is back-illuminated 
(BI: XIS1) (Koyama et al. 2007). The XIS was operated in 
the normal clocking mode (8 s exposure per frame), with 
the standard 5x5 and 3x3 editing mode. We processed 
the XIS data using the "xispi" and "makepi" ftools task 
and the CALDB files of version 2009-08-13. The XIS data 
were then cleaned by assuming thresholds on the Earth el- 
evation angle of > 5° and the Day-earth elevation angle 
of > 20°. We also discarded data for time since passage 
through the South Atlantic Anomaly of less than 436 sec. 
We created a light curve for each sensor over 0.5-2 keV 
binned to 540 sec but there are no anomalous event rate 
greater or less than 3er around the mean. After the above 
screening, the remaining exposure was 102.3 ksec for both 
FIs and BI. Event screening on cut-off rigidity was not 
performed. 



3. Spectral Analysis with Chandra and Suzaku 
Data 

3.1. Definition of regions for Spectral Analysis 

The left and right hand panels of figure 1 show the X- 
ray source images from Chandra and Suzaku overlaid with 
optical contours constructed from Digitized Sky Survey 
data. Because the starburst region rises about 1.3 kpc 
(~15") from the center (Cecil et al. 2002), we extracted 
spectra from the following regions of Chandra and Suzaku 
images. For Chandra, there are three regions centered on 
the nucleus at coordinates (R.A., Dec.) = (10 h 01 m 58f51, 
+55°40'49"4). 1) a circle with a 0.5' radius (2.25 kpc; 
hereafter '0.5' circle'); 2) an inner ring between 0.5'-l' 
(2.25-4.5 kpc; '0.5'-l' ring'), and 3) an outer ring between 
l'-2' (4.5-9 kpc; 'l'-2' ring'). The background region for 
the Chandra analysis is a ring with a radius of between 
2'-3' as shown in figure 1 (left). The count rate in this 
background ring is less than 5% of that measured within 
the central 2'. 

For Suzaku, the source region is a circle with a 
4' radius ('4' circle') centered on coordinates (R.A., 
Dcc.) = (10 h 01 m 59?69, +55°39'59"8). This is much larger 
than for Chandra owing to Suzaku's half-power diame- 
ter of ^2'. The background region is the entire XIS field 
of view, excluding the source region and two additional 
circular apertures with radii of 2' each, centered on back- 
ground field objects as shown in figure 1 (right). Spectral 
fitting was performed with XSPEC 12.6 (Arnaud 1996). 

3.2. Point Sources with Chandra 

The low angular resolution of Suzaku, ~2', prevents 
separation of point sources (PSs) but Chandra is ideal for 
this purpose. The PSs were detected by the "wavdetect" 
CIAO tool task, with parameter 'sigthresh' set to its de- 
fault value of 10~ 6 in the 2-10 keV energy range. We accu- 
mulated the spectra of all detected PSs within 4' of the nu- 
cleus using regions determined by "wavdetect" , and fitted 
their sum with the following model: phabspg x (power- 
lawpg + Gaussian). Here, phabspg represents photoelec- 
tric absorption. The integrated spectra from PSs are usu- 
ally approximated by a power-law model, denoted above 
by power-lawps (e.g., Xu et al. 2005). Since Iyomoto et 
al. (2001) detected a strong Fe-K line at 6.4 keV, we also 
include a Gaussian above. This model reproduces the ac- 
cumulated spectrum well with x 2 /d.o.f.=46/36. The de- 
rived parameters are listed in table 1 and the integrated 
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Table 2. Best-fit parameters for the background regions of 
Suzaku with apec components + PL model.* 



cm 



Parameters 

IVk (10 
Tcxb 

fcTMWH (keV) 
Abundance (solar) 
fcT LHB (keV) 
Abundance (solar) 
Norm hHB /Norm MWH 

X 2 /d.o.f. 



0.795 (fix) 
1 4fi+ - 08 

0.30±0.03 

1 (fix) 
0.10±0.01 

1 (fix) 

1 9 oc+8.30 

333/303 



* The apec components for spectra in the background region of 
NGC 3079 with absorbed MWH and LHB components for 
the Galactic emission, and a PL model for the CXB. 

spectrum is shown in figure 2. The index of power-lawps 
is 1.92+o 28> a typical value for AGN spectra. The central 
line energy of Gaussian model is 6.30+o;o3 keV, which is 
consistent with that of Iyomoto et al. (2001) and Fukazawa 
ct al. (2011). 

3.3. Estimation of Background Spectra with Suzaku 

For the Suzaku data analysis, we performed simulta- 
neous fitting of the source and background spectra, in- 
stead of subtracting the background as with Chandra. 
The reason is because of the extended tail of point 
spread function and the larger vignetting in the case of 
Suazaku, as compared to Chandra (Scrlemitsos et al. 
2007). Therefore, we first fitted the spectra extracted 
from the Suzaku background region. Since the NXB is 
already subtracted, the background components consist 
of Cosmic X-ray Background (CXB) and Galactic ther- 
mal emission. We assumed a power-law model for the 
CXB component, and a two-temperature thermal plasma 
model for Galactic emission. One of the plasmas is the 
Milky Way Halo (MWH) and the other is the sum of 
the solar wind charge exchange (SWCX) emission and the 
Local Hot Bubble (LHB) (Yoshino et al. 2009). 

The fitting was performed with the following model; 
phabs G x (power- lawcxB + apecMWri) + apecLHB- 
phabsg represents photoelectric absorption with column 
density fixed to the Galactic value 0.795 x 10 20 cm -2 in the 
direction of NGC 3079. The apccMWH and apecLHB com- 
ponents represent thin thermal plasmas using the APEC 
model (Smith et al. 2001) from MWH and LHB (and 
SWCX), respectively. The redshift and metal abundance 
are fixed at zero and one solar, respectively. We simulta- 
neously fitted the spectra from BI and FI CCDs in the 0.4- 
5.0 and 0.5-5.0 keV range, respectively. The derived pa- 
rameters are listed in table 2. This model reproduced the 
spectra approximately, with x 2 /d.o.f.=333/303. Because 
the photon index and flux of power-lawcxB and temper- 
atures of apec models are consistent with typical values 
(e.g. Kushino et al. 2002; McCammon et al. 2002; Yoshino 
et al. 2009), these models represent the background spec- 
tra well. We adopt these parameters for our background 



model. 

3.4- Simultaneous Fit of Chandra and Suzaku Data 

We fitted simultaneously the spectra of the three 
regions (0.5' circle, 0.5'-l' ring, and l'-2' ring) of 
Chandra background-subtracted data, along with the 
spectra from the 4' circle and background region of 
Suzaku. Each Chandra spectrum was fitted with 
phabsc x vapeco.5' or o.5'-i' or i'-2'- The phabsc fac- 
tor is fixed to the Galactic value 0.795xl0 20 cm~ 2 . 
The vapeco.5' or o.5'-i' or i'-2' model represents thermal 
plasma emission in the ISM of NGC 3079, corresponding 
to the 0.5' circle, 0.5'-l' ring, and l'-2' ring regions of 
Chandra, respectively. We divided the metals into the 
groups of O, Ne, (Mg & Al), (Si, S, Ar, & Ca) and (Fe 
& Ni) based on the metal synthesis mechanism of SN, al- 
lowing each group to vary. However, the value of the Si 
group is fixed to 1 solar beyond the 1' circle region. 

For the Suzaku 4' circle spectra, we included all three 
Chandra components in the models. We fitted the spec- 
tra of Suzaku 4' circle with following model: phabsc 
x (vapeco.5' + vapeco.5'_i' + vapcci'_2') + phabsps x 
(power-lawpg + Gaussian) + background models. The 
first term represents the three components of Chandra re- 
gions, which were linked to that of corresponding each 
Chandra region. The accumulated emission from PSs in- 
cluding the nucleus is modeled as phabsps (power-lawps 
+ Gaussian). The phabsps and index of power-lawps are 
fixed to the best fit values derived in subsection 3.2. The 
parameters of the background models were linked between 
the 4' circle and Suzaku background regions. An energy 
range of 0.4-4.0 keV is used for the 0.5' circle, 0.5'-l' ring, 
and l'-2' ring of Chandra, and from 0.4 (BI) or 0.5 (FI) 
to 5.0 keV and 6.0 keV for the background and 4' circle of 
Suzaku, respectively. Because there are known calibration 
issues of the Si edges in all XIS CCDs, the energy range of 
1.84-1.86 keV is ignored (Koyama et al. 2007) in Suzaku 
spectra. 

The fitted spectra are shown in figure 3 and the derived 
parameters are listed in table 3. This model reproduced 
these spectra approximately, with x 2 /d.o.f.=885/805. 
The derived temperatures of thermal models are 0.65, 
0.45, and 0.24 keV in 0.5' circle, 0.5'-l' ring, and l'-2' ring 
regions, respectively. Several emission lines seen around 
0.5-0.6 keV, 0.6-0.7 keV, -1.3 keV, and -1.8 keV are 
identified with Ka lines of O VII, O VIII, Mg XI, and Si XIII, 
respectively The emission bump around 0.7-1 keV cor- 
responds to the Fe-L complex, as well as to K-lines from 
Ne IX and Ne X. The metal abundances in all three regions 
have very similar best-fit values, though with large uncer- 
tainties. In particular, the lowest temperature (0.24 keV) 
component makes the faintest contribution to the overall 
spectrum, so the low O abundance in this region should 
have a large associated systematic error. Meanwhile, the 
O abundance derived from highest temperature (0.65 keV) 
bright component is more reliable. 

In order to examine the abundance ratios rather than 
their absolute vales, we calculated the confidence contours 
between the abundance of various metals (O, Ne, Mg, and 
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Table 3. Best-fit parameters for the 0.5' circle, 0.5'-l' ring, l'-2' ring, and 1' circle regions of Chandra and 4' circle of Suzaku 
with vapec components, PL model, and Gaussian model. 



Parameters 



Chandra Suzaku Chandra Suzaku 

1' abundance linked 







0.5' circle 








0.5' circle 






(10 20 cm- 2 ) 


0.795 (fix) 


0.795 


(fix) 


0.795 (fix) 


0.795 (fix) 


kT 


(keV) 


0-65™ 


linked 


to 


lett 


0.63±0.04 


linked to left 





(solar) 


0.6118$ 


linked 


to 


leit 


n 5n +0 - 24 

u - ou -0.19 


linked to left 


Ne 


(solar) 


1 oo +0 - 46 

-L.UU_ .3 4 


linked 


to 


left 


n qq+0-27 

U.»»_ Q 23 


linked to left 


Mg, AI 


(solar) 


n 4fi+ - 28 

u -^°-0.22 


linked 


to 


left 


0.48t°' 2 ? 


linked to left 


Si, S, Ar, Ca 


(solar) 


S7+ 29 

u -°' -0.22 


linked 


to 


left 


42+ 24 

u -^ z -0.20 


linked to left 


Fe, Ni 


(solar) 


24+ - 08 

u - z ^-0.05 


linked 


to 


left 


n 99+O.O5 
u - zz -0.04 


linked to left 


Norm* 


(io- 4 ) 


1 qn+o.38 

1.3U_o 39 


linked 


to 


left 


1 q2+ - 34 

i-»^_0.31 


linked to left 


O/Fc 


(solar) 


2 56+ 179 

z - OD -1.36 


linked 


to 


left 


9 1 O + 1.03 

z - lo -0.77 


linked to left 


Ne/Fe 


(solar) 




linked 


to 


left 


4 50+ 1 - 31 

^■ ou -1.00 


linked to left 


Mg/Fe 


(solar) 


9 i 7+1.53 
Zl1 ' -1 41 


linked 


to 


left 


9 i 7+0.91 
' -0.77 


linked to left 


Si/Fe 


(solar) 


1 64+ 074 


linked 


to 


left 


JU -0.94 

0.5'-l' ring 


linked to left 


total counts 


(count) 


2327 
0.5'-l' ring 










kT 


(keV) 


4^+ 07 


linked 


to 


left 


0.4«_o. 09 


linked to left 





(solar) 


41 +°- 41 


linked 


to 


left 


linked to 0.5' 


linked to left 


Ne 


(solar) 


1 n7+°- 60 

1 - u '-0.37 


linked 


to 


left 


linked to 0.5' 


linked to left 


Mg, AI 


(solar) 


0.67+°-^ 

— u.oy 


linked 


to 


left 


linked to 0.5' 


linked to left 


Si, S, Ar, Ca 


(solar) 


1 (fix) 


linked 


to 


left 


linked to 0.5' 


linked to left 


Fe, Ni 


(solar) 


23+ 10 

u.zo_o 07 


linked 


to 


left 


linked to 0.5' 


linked to left 


Norm* 


(io- 4 ) 


1.11±0.32 


linked 


to 


left 


1 14+0.26 
i - i ^-0.22 


linked to left 


O/Fc 


(solar) 


i fi n +2 - 40 

A - ou -0.80 


linked 


to 


left 


linked to 0.5' 


linked to left 


Ne/Fe 


(solar) 


4 60+ 216 
^• DU -l.00 


linked 


to 


left 


linked to 0.5' 


linked to left 


Mg/Fe 


(solar) 


q i q+1.87 
°- 10 -1.84 


linked 


to 


left 


linked to 0.5' 


linked to left 


total counts 


(count) 


3898 
l'-2' ring 








l'-2' ring 




kT 


(keV) 


24+ - 03 

u - z ^-0.02 


linked 


to 


left 


0.24±0.03 


linked to left 





(solar) 


0-30t°; 27 


linked 


to 


left 


o.3ii_;_| 


linked to left 


Ne 


(solar) 


L -°'-0.69 


linked 


to 


left 


i 7n +1 - 56 

1 -' v -0.71 


linked to left 


Mg, AI 


(solar) 


Q7+4.56 
z - 8 '-1.93 


linked 


to 


left 


o 90+3.34 
°- zo -2.07 


linked to left 


Si, S, Ar, Ca 


(solar) 


1 (fix) 


linked 


to 


left 


1 (fix) 


linked to left 


Fe, Ni 


(solar) 


i 99+1.44 
1 - zz -0.59 


linked 


to 


left 


1 - zu -0.62 


linked to left 


Norm* 


(io- 4 ) 


i-6o±S:g 


linked 


to 


left 


1 ^q+0-62 

1 - oy -0.64 


linked to left 


O/Fe 


(solar) 


0.24t°; 4 I 


linked 


to 


left 


24+ 29 

u - z ^-0.13 


linked to left 


Ne/Fe 


(solar) 


l-OU_o 40 


linked 


to 


left 


J - aD -0.61 


linked to left 


Mg/Fe 


(solar) 


QK+1-46 
z -°°-1.40 


linked 


to 


left 


2A2±\f 7 


linked to left 


^H PS 


(10 20 cm" 2 ) 




115 (fix) 




115 (fix) 


Tps 






1.92 (fi: 


<0 




1.92 (fix) 


total counts^ 


(count) 


13291 


63313 






X 2 /d.oi. 






885/805 




886/809 



* Normalization of the vapec component divided by the solid angle, Q v of the source image in "xissimarfgen" , Norm = J n e nndV / (4ir (1 + 



z) 2 Dl)/n v xio- 



wherc n c and njj is the electron and proton density, respectively, z is the redshift of the source, and 



Dj± is the angular diameter distance to the source. 

' Total counts are not subtracted background emission and is accumulated of all XIS detectors in Suzaku. 
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Fig. 3. Simultaneous fitting of Chandra and Suzaku data. Top panels: Background-subtracted ACIS-3 spectra of the 0.5' 
circle (left), 0.5'-l' ring (middle), and l'-2' ring (right) region, shown without removal of instrumental response. Magenta, 
orange, and blue lines denote the ISM component (vapec) each region. Bottom panels: NXB-subtracted XISO (black) and 
XIS1 (red) spectra of the 4' circle of NGC 3079 (left panel), and those of the background region (right panel), shown with- 
out removal of instrumental response. Black and red lines show the best-fit model for XISO and XIS1, respectively. For 
simplicity, only the model components for XIS1 spectra are shown. Green line shows accumulated emission from point 
sources (absorbed power-lawps). The magenta, orange, and blue lines are linked to Chandra spectra. The gray lines are 
the Galactic background emission (apec^wH an d a P ec LHB) ancl the CXB components, respectively. The background com- 
ponents are common between the on-source and background spectra, but scaled to the respective data accumulation area. 



Si) to that of Fe, in order to derive abundance ratios. The 
results are shown in table 3 and figure 4, where we also 
indicate 90%-confidence abundance of these metals rela- 
tive to Fe with black dashed line. The elongated shape 
of the confidence contours indicates that the relative val- 
ues were determined more accurately than the absolute 
values. Because the abundance patterns between the 0.5' 
circle and 0.5'-l' ring are consistent with each other, we 
fitted the Chandra and Suzaku spectra again, this time 
setting the abundances to be common between these two 
regions. The derived parameters agree with the results in 
previous free fit, as summarized in table 3 and figure 4, but 
now the abundances in the inner regions arc better con- 
strained. In figure 4, we only plotted the 90%-confidence 
abundance of these metals relative to Fe with red dashed 
line. The calculated abundance ratios of O/Fe, Ne/Fe, 
Mg/Fe, and Si/Fe, are 2.13^;^, 4.50±J;gJ, 2.17^;^, and 
1.90i ' g 4 within the 1' circle, respectively. 

3.5. Systematic Uncertainties in the Abundance Ratios 

In Chandra analysis, we extracted the spectra after ex- 
cising the PSs. We fitted the spectra by setting the power- 
law model normalization within the 0.5' circle region to be 
free, in order to investigate the effect of the emission from 
PSs below the detection threshold. The flux upper limit 
of the power-law model corresponds to <20% of the total 



flux in the central region. Furthermore, the parameters of 
vapec model do not change significantly. 

In order to investigate systematic uncertainties for the 
abundance ratios, we set free the absorption parameter, 
phabsc, because additional absorption may arise along 
the edge-on line-of-sight within the ISM of NGC 3079. We 
tried setting free the absorption affecting only the central 
0.5' circle region, where intrinsic absorption should be the 
highest. The resultant absorption was constrained to an 
upper limit of 8.4 xlO 20 cm~ 2 . None of the other param- 
eters changed significantly, implying that the systematic 
effect of the absorption uncertainty is almost negligible. 

Next, we investigated abundance trends amongst the 
a-elcmcnts. In particular, we set the abundances of Ne 
and Mg to be equal to that of O in each region because 
we expect similar distributions among these. In the 0.5' 
circle and 0.5'-l' ring regions, the ratios of O/Fe, Ne/Fe, 
and Mg/Fe are consistent to those of each region, as in 
subsection 3.4 above. On the other hand, in the l'-2' ring 
region, the ratio of O/Fe increased to near solar. As men- 
tioned in subsection 3.4, the O abundance of the lowest 
temperature component may have a large systematic error 
owing to its faintness. 

The Fe-L lines have systematic uncertainties related to 
unknown atomic physics, as discussed in Matsushita et al. 
(2000). When convolving with the XIS response matrix, 
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Fig. 4. Confidence contours between metal (O, Ne, Mg, and Si) and Fe abundances, determined for 0.5' circle re- 
gion. The black, red, and green contours represent 68%, 90%, and 99% confidence ranges, respectively. The black 
dashed lines indicate 90%-confidcnce abundances of these metals relative to Fe. The red dashed contours show the 
90% confidence ranges for when the abundances in the 0.5' circle and 0.5'— 1' ring regions are linked to each other. 



differences in the Fe-L lines between APEC (Smith et al. 
2001) and MEKAL (Mewe et al. 1985) models are typi- 
cally of order 10-20% (Matsushita et al. 2000). Therefore, 
we re-fitted the spectra in the same way as above, but 
with an inclusion of 15% systematic errors over the en- 
ergy range of 0.75-1.2 keV. The x 2 /d.o.f. improves to 
603/623, with null-hypothesis probability of ~0.7. The 
derived parameters in this case are completely consistent 
to those of subsection 3.4. Therefore, these results are use- 
ful to assess whether the temperatures and metal abun- 
dances in the ISM were reasonably determined or not. 
Furthermore, the APEC code has recently been updated 
to version 2. The ratios among L-lincs complex of Fe have 
changed significantly from version 1. We have fitted the 
spectra with version 2 of APEC code, and found that all 
resultant parameters are consistent with those of table 3, 
with x 2 /d.o.f.=893/805. In this statistic, there are no dif- 
ferences between results in version 1 and 2 of APEC code. 



4. Discussion 

4--1- Spatial Metal Abundance Patterns 

We have successfully derived the spatial metal abun- 
dance patterns in the hot ISM of NGC 3079, using both 
Chandra and Suzaku for the first time. The combination 
of these missions enables us to optimally utilize the advan- 
tages of both good spatial and good spectral resolution at 
high sensitivity. Figure 5 shows the abundance patterns 
in the 0.5' circle, 0.5'-l' ring, l'-2' ring, and 1' circle re- 
gions. The abundance patterns of NGC 3079 can now be 
compared with the solar abundance table (Lodders 2003) 
and those expected from SN II and SN la yields. The SN 
II yields computed by Nomoto et al. (2006) are used here. 
These refer to an average over the Salpeter initial mass 
function for stellar masses from 10 to 50 Mq, with a pro- 
genitor metallicity of Z = 0.02. The SN la yields from the 
W7 model of Iwamoto et al. (1999) were adopted. Within 
a 4.5 kpc region (i.e. a circle with radius of 1'), the abun- 
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Fig. 5. Abundance ratios of O, Ne, Mg, and Si to Fe for the best- fit model of NGC 3079 derived by simultaneously 
fitting the Chandra and Suzaku spectra. Magenta, orange, blue, and purple data show the abundance patterns of 0.5', 
0.5'— 1', 1'— 2', and 1', respectively. Solid, dot-dashed, and dotted lines represent the number ratios of metals to Fe for so- 
lar abundance, for SN II, and for SN la products (Lodders 2003; Iwamoto ct al. 1999; Nomoto et al. 2006), respectively. 



Table 4. Spatial scales of well-known starburst galaxies, in which abundance patterns are SN II— like. The distance scale is measured 
perpendicular to the disk in each case. 



Galaxy 


NGC 3079 NGC 253* 


NGC 4631 T 


M 82 


D 


(kpc) 4.5 9 


13 


11* 



* Bauer et al. (2008) t Yamasaki et al. (2009) * Tsuru et al. (2007) 
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dance pattern is consistent to that of SN II within the 
uncertainties. We thus infer this region to be significantly 
enriched by SN II metal yields, undoubtedly related to 
the starburst activity in this system. On the other hand, 
the abundance pattern is not SN II— like beyond 4.5 kpc, 
presumably because the starburst has had little impact 
here. This emission comes from the intrinsic galactic halo 
and is not enriched by the starburst outflow yet. 

Using AKARI, Yamagishi et al. (2010) detected ex- 
tended emission from hot dust (~30 K) in the center re- 
gions (4 kpc) of NGC 3079, possibly associated with the 
starburst activity. The emission ratios of ionized to neu- 
tral polycyclic aromatic hydrocarbons are higher in the 
center (^2 kpc) than in the disk region. Both facts indi- 
cated that star formation is relatively active in the center 
region, consistent with our results based on X-ray data. 

4-2. Time Scale of Starburst Activity 

The total infrared luminosity of NGC 3079 indicates a 
rather low star-formation rate, and measured gas-to-dust 
mass ratios over the galaxy are unusually high, which led 
Yamagishi et al. (2010) to suggest that NGC 3079 is still 
in a rather early phase of the starburst, with nuclear star 
formation activity expected to increase in the near future. 
Derived abundance patterns indicate that the ISM 
within at least 4.5 kpc from the center contains metals 
synthesized by SN II. Using this fact, the approximate 
time (t) elapsed since starburst activity began can be 
estimated by computing the distance (R) travelled by the 
wind according to the numerical hydrodynamic starburst 
wind simulations of Tomisaka & Bregman (1993). R can 
be written as 



R(t) ~ 31 kpc ( 



Uhalo 



-1/5 



(: 



2 x 10- 2 cm- 2 ' 



\3/5 



*3 x 10 42 erg s -1 ' '50 Myr ; 

where L is the luminosity and rihaio is the particle den- 
sity of hot plasma (Tomisaka & Bregman 1993). From 
the Suzaku spectra, we measure a 0.3-5 keV total un- 
absorbed ISM luminosity of 7.33 xlO 39 erg s _1 . The 
plasma (electron) density is derived from the normaliza- 
tion of vapec model, Norm = J n c nndV / (4tt (1 + z) 2 D\) 
xlO -14 cm -5 arcmin -2 , where £>a is the angular distance 
to the source. The normalization of vapec model within 
1' circle region is 3.0 x 10~ 4 enr 5 arcmin -2 (table 3). 
Considering the ratio of proton to He in plasma for SN 
II— like abundances, we use a particle density of 5 x 10~ 3 
cm -3 , as obtained from the emission measure defined by 
/ 2.67n e ri}idV. With these assumptions, we derived the 
time for the wind to travel a distance of 4.5 kpc as 10 Myr. 
In combination with the 4 kpc extent of the infrared- 
emitting starburst region, 10 Myr may be considered as 
an upper limit for the age of the starburst. 

In two-dimensional numerical hydrodynamical simula- 
tions incorporating more realistic spatial plasma density 
distribution (e.g. Tomisaka & Bregman 1993, Strickland 
& Stevens 2000), starburst-driven galactic winds can reach 



~ 5 kpc after 5 Myr < t < 10 Myr. Our calculation is 
comparatively crude (assuming spherical symmetry and a 
uniform density), but the results are qualitatively consis- 
tent. 

4-3. Comparison of the Starburst Phases among 
Starburst Galaxies 

In a few other well-known starburst galaxies, 
NGC 4631, M 82, and NGC253, the metal abundances 
of the outflow regions also show a-elemcnt enhance- 
ments relative to Fe, consistent with enrichment by SN 
II (Yamasaki et al. 2009, Tsuru et al. 2007, Konami et al. 
2011, Mitsuishi 2012). The central region of NGC 3079 
shows a similar pattern. In table 4, we summarize the 
spatial extents of the regions with SN II— like patterns in 
these galaxies, measured perpendicular to the disk in each 
case. Spatial scales of order ~10 kpc in NGC 4631, M 82, 
and NGC253 indicate that more than 10 Myr have passed 
since starburst activity began, according to the simula- 
tions of Tomisaka & Bregman (1993). By comparison, 
the smaller spatial extent and age upper-limit inferred in 
subsection 4.2 indicate that NGC 3079 is in an earlier 
phase of activity. Though the SN II— like abundance pat- 
terns have been observed through starburst galaxies, the 
different spatial distributions of those patterns correspond 
to their starburst activity phases. 

Though the temperature of dust in the center of 
NGC 3079 is about ~ 30 K, in the central (<1 kpc) region 
of M 82 at least two components with higher dust temper- 
atures of ~45 K and 160 K have been found (Yamagishi 
et al. 2010, Gandhi et al. 2011). These may result from 
prolonged heating associated with the longer period of 
starburst activity in M 82, as compared to NGC 3079. 
Copious production of dust in more mature starbursts 
like M 82 also ought to absorb a larger fraction of the 
starburst radiative output and thus become more heated. 
NGC 3079 instead displays a comparatively low dust-to- 
mass ratio at present. For NGC 253, Engelbracht et al. 
(1998) concluded that the source is in a late starburst 
phase, based upon detection of extended infrared emis- 
sion. Modeling of this emission implied an age of 20 to 
30 Myr. In NGC 4631, the ratios of 12 CO(2-l)/ 12 CO(l- 
0) are similar to the average value seen in the nuclei of 
normal galaxies (Golla & Wielebinski 1994). The authors 
used this fact to suggest that source was undergoing mild 
central star formation. But these ratios could also be con- 
sistent with a late phase of starburst activity. 

The abundance pattern beyond 4.5 kpc in NGC 3079 
is rather consistent with solar (Lodders 2003). The disk 
region in the starburst galaxy NGC 4631 similarly shows 
solar abundance. It is interesting that the hot ISM in 
NGC 4258, a spiral galaxy without starburst activity 
(Konami et al. 2009), also behaves similarly. These facts 
suggest that the metallicity of the ISM both before, and 
after an era of intense starburst activity, may look quite 
similar. 
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5. Conclusion 

We have performed X-ray spectral analysis of NGC 3079 
with Chandra and Suzaku. The temperatures in the ISM 
thermal plasma are found to show a negative radial gradi- 
ent. In the central (<4.5 kpc) region, the metal abundance 
patters are consistent with those synthesized by SN II. On 
the other hand, abundances in the outer region (>4.5 kpc) 
are closer to solar. These indicate that the metal abun- 
dances in the center have been enriched with SN II yields, 
but that the outer region has not yet been influenced by 
the starburst. This is consistent with starburst activity in 
NGC 3079 being in a relative-early evolutionary phase. 

We thank the referee for providing valuable comments. 
We gratefully acknowledge all members of the Suzaku 
hardware and software teams and the Science Working 
Group. SK is supported by JSPS Research Fellowship for 
Young Scientists. 
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